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ABSTRACT : Leptin is a 16 kD protein, synthesized by adipose tissue and is involved in regulation 
of feed intake, energy balance, fertility and immune functions. Present study was undertaken with 
the objectives of sequence characterization and studying the nucleotide variation in leptin gene in 
Murrah buffalo. The leptin gene consists of three exons and two introns which spans about 18.9kb, of 
which the first exon is not transcribed into protein. In buffaloes, the leptin gene is located on 
chromosome eight and maps to BBU 8q32. The leptin gene was amplified by PCR using oligonucleotide 
primers to obtain 289 bp fragment comprising of exon 2 and 405 bp fragment containing exon 3 of 
leptin gene. The amplicons were sequenced to identify variation at nucleotide level. Sequence 
comparison of buffalo with cattle reveals variation at five nucleotide sequences at positions 983, 
1083, 1147, 1152, 1221 and all the SNPs are synonymous resulting no in change in amino acids. Three 
of these eight nucleotide variations have been reported for the first time in buffalo. The results indicate 
conservation of DNA sequence between cattle and buffalo. Nucleotide sequence variations observed 
at leptin gene between Bubalus bubalis and Bos taurus species revealed 97% nucleotide identity. 
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INTRODUCTION 
Leptin, a 16 kD a protein that is synthesized 


leptin on milk production. Several single 
nucleotide polymorphisms (SNP) have been 


by adipose tissue is involved in regulation of 
feed intake, energy balance, fertility and 
immune functions (Fruhbeck et al. 1998). Milk 
production trait is a quantitative trait and 
polygenic in inheritance. Since milk production 
trait is directly related to feed intake and energy 
balance, it is obvious that there is effect of 


reported in the coding region of the gene 
(Konfortov et al. 1999, Orrt et al. 2007, 
Liefers et al. 2002) found that some SNP are 
associated to production, live weight, energy 
balance, feed intake and fertility in Holstein 
heifers. Leptin is a 167-amino acid protein 
produced by the leptin gene (leptos means thin). 
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Leptin is also responsible for the regulation of 
body weight and energy homeostasis and is one 
of the potential biomarker for high performance 
leading to better adaptability and productivity. 
(Friedman and Halaas 1998). The Leptin gene 
consists of three exons and two introns which 
spans about 18.9 kb, of which the first exon is 
not transcribed into protein. Leptin gene has 
been mapped to chromosome 7 in human 
(Green et al. 1995) chromosome 4 (BTA4) in 
bovine (Stone et al. 1996) and chromosome 8 
(BBU 8q32) (Vallinato et al. 2004) in buffalo. 
Several SNPs have been previously identified 
in introns and exons of leptin among different 
breeds of cattle. 

Leptin plays a critical role in regulating and 
coordinating energy metabolism (Friedman and 
Halaas 1998). It regulates the metabolism of 
key tissues involved in the storage and 
dissipation of energy (Banks et al. 2000). 
Therefore, leptin may be important in 
regulating metabolic adaptation of nutrient 
partitioning during the energy-consuming 
processes of pregnancy and lactation (Moschos 
et al. 2002). Since buffalo contributes about 
54 per cent of the total milk produced in India. 
Although the economic importance of 
buffaloes has always been known, yet very little 
work has been carried out to exploit the 
genetic potentials of this animal. Hence the 
present study was undertaken with the objective 
to identify sequence variation in leptin gene of 
buffalo. 


MATERIALS AND METHODS 

Sample Collection and DNA preparation 

Blood samples were collected in vacutainers 
(Bacton-Dickinson vacutainer system) 
containing sodium EDTA as an anticoagulant 
from Murrah buffaloes maintained at National 
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Dairy Research Institute, Karnal. Genomic 
DNA was isolated from blood samples 
following phenol-chloroform extraction 
method described by Sambrook and Russel 
(2001) with slight modifications. DNA was 
dissolved in TE buffer at 60°C in water bath 
for 2 hrs to dissolve pellet properly in buffer 
and was kept in refrigerator for further use. 
Quality of DNA was checked through 
spectrophotometry. DNA samples with O.D. 
ratio between 1.7 and 1.9 were considered as 
good and used for further study. The samples 
beyond this range were re-extracted by Phenol- 
chloroform extraction method. DNA quality 
was also checked by running the sample in 0.8 
percent agarose gel electrophoresis. The quality 
of DNA was checked on 0.8% agarose and 
quantity by UV spectrophotometer at A,./A,.., 
nm. 


DNA Amplification 

The PCR amplification using primers 
(Adikari 2006) for exon 2 and 3 of leptin gene 
was carried out in a programmable thermal 
cycler (PTC 200, MJ Research) using the 
following program: For the amplification of 2" 
exon, 25 ul of PCR reaction mix contained 
3.000 ul genomic DNA (50 ng/ul), 0.600 ul 
primers each (100 pM/ul), 12.5 ul Fermentas 
Master Mix™ (2X) and 8.3 ul double distilled 
H,O. , Initial denaturation at 93°C for | min, 
denaturation at 93°C for 1 min, annealing at 
59°C for 30 sec, extension at 72°C for 1.3 min 
were carried out for 36 cycles. 

For the amplification of 3rd exon, 25 ul of 
PCR reaction mix contained 3.000 pl genomic 
DNA (50 ng/ul), 0.600 wl primers each (100 
pM/ul), 12.5 ul Fermentas Master Mix™ (2X) 
and 8.3 wl double distilled H,O, Initial 
denaturation at 93°C for 1 min, denaturation 
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at 93°C for 1 min, annealing at 56°C for 30 
sec, extension at 72°C for 1.3 min were carried 
out for 36 cycles. The PCR products were 
loaded on 1.5% agarose to confirm the 
amplification of target region using 100 bp 
ladder as a marker and were visualized by 
ethidium bromide staining, which was added 
to the agarose gel @ 1 ul /100 ml of gel. The 
agarose gel electrophoresis was performed in 
1X buffer at 100 volts for 45 minutes and PCR 
products were visualized on UV 
transilluminator and photographed with gel 
documentation system. The primers designed, 
regions amplified, annealing temperatures and 
product sizes are given in Table]. 


Nucleotide Sequencing 
The PCR products were concentrated to 50 
ng/u by pooling several tubes to precipitate by 


of the DNA to UV light). The minimum agarose 
slice was transferred to a 1.5 ml microcentrifuge 
or screw cap tube and then purified by using 
commercially available gel extraction kits 
(Qiagen). Quantification was done by loading 
one ul of eluted sample in 1% Agarose gel and 
comparing with standard molecular marker (Phi 
X 174 DNA ladder or 100 bp DNA ladder). Only 
samples with good concentration (>50 ng/l) 
were selected and subjected to sequencing. 
Amplified PCR products from each set of 
primers were subjected to custom DNA 
sequencing from both ends (5' and 3' ends). from 
M/s. Chromous biotech, INDIA. 


Sequence Data Analysis 

Sequence data were analyzed by using 
Chromas (Ver. 1.45, http://www.technelysium. 
com.au/chromas.html). Sequence data from 


Table 1:Description of primers used and the amplified product of different loci studied 


Primer sequence 5’-3’ Tm Amplified product length 


F-5’- GGT GGT AAC GGA TCA CAT GG -3’ 
R-3'- CCA CGG TTC TAC CTC GTC TC 3’ 


F-5’- GCA TAG CAG TCC GTC TCC TC -3’ 
R-3’- TTC CCT GGA CTT TGG GAA G -3’ 


the isopropanol procedure. In order to obtain 
clean fragment for sequencing, the PCR 
products were separated by electrophoresis in a 
TAE agarose gel containing ethidium bromide 
using standard protocols. The desired PCR 
product band was excised using a clean, sterile 
razor blade or scalpel (band was visualized in a 
medium or long wavelength (e.g. 300 nm) UV 
light, and excised quickly to minimize exposure 


mi ad 289 bp fragment 


contaning exon II 


56°C 


405 bp fragment 


containing exon II 


variants of different regions were subjected to 
multiple alignments (DNASTAR, Clustal W) 
for identifying the SNPs. 


Database Search 

The database search of sequences for a 
possible match to the DNA sequence of leptin 
gene was conducted using the BLAST 
algorithm available at the National Center for 
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Biotechnology Information (NCBI). Translated 12 new single nucleotide polymorphism (SNP); 
protein sequences of leptin genes in different two polymorphic sites G3333A and G3441A 
species namely Bos taurus, Bos indicus, Ovis previously identified in exon 3 by Vallinoto et 
aries, Capra hircus and Homo sapiens were al. (2004) and originally referred to as 276-SNP, 


also subjected to BLAST algorithm. 384-SNP and one of the new SNP, C1221T 
occur in exon 2 and three G3195A, G3318A 
RESULTS AND DISCUSSION and G3434A in exon 3 while SNPs T1015C, 


Sequence data were analysed using chromas C1071T, G1072A, T1081C, T1131G, T1143C, 
(Ver.1.45, http://www.technelsium.com.au / — T1145G and A1151G are all in intronic region. 
chromas.html). ClustalW multiple alignments —_ They also reported a mutation at G3434A in the 
with Bubalus bubalis sequence reference exon 3, resulting in Arg to Gln substitution in 
sequence NCBI GI 46254758) revealed the codon 159 but it has been observed only 
nucleotide changes at 5 positions in the gene. one Italian Buffalo while all other animals were 
These variations are found in nucleotide | homozygote GG genotype. In exon 3, the SNP 


Table 2: SNPs Identified in Leptin Gene (Murrah buffaloes) 


Position in Nucleotide SINP type Amino 
aligned change add 
data change 


sequences at positions G983A in intron 1 which = G3318A was found polymorphic only in the 
is heterozygous; at C1083T-homozygous , Egyptian Buffaloes, while it resulted 
T1147G- heterozygous, Al1152G-homozygous monomorphic in Italian Buffaloes. Position 
and C1221T- heterozygous mutations in exon G3441A was found monomorphic in both 
2; atT1371C in intron 2 homozygous mutation; | Egyptian and Italian Buffalo. 
(Fig.I). However, all these nucleotide changes Fitzsimmones et al. (1998) reported 4 SNPs 
are synonymous i.e. there isnochangeinamino —‘ in exon 2 in Bos taurus. Whereas, Haegman et 
acids. (Table 1) Nucleotide sequence data and al. (2002) found 2 SNPs in exon 2 in Bos 
annotation details are given in Fig.2 . taurus. Adikari (2006) has also reported 
Orru et al. (2007) while investigating River nucleotide changes at position 11, 365, 369 and 
Buffalo and two Egyptian River Buffalo reported 371 in intron 2 region. Six SNPs in exon 2 of 
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Fig.1: Clustal W Alignment with Reference sequence of Leptin Gene (Exon-2) in Murrah Buffalo 


GGATTAGTGCCTTICATTACTGICATTTCTAGACAATGAATIGICTITGAGGAGATGATAGCCATGGCAGACAGC 
) TIGTTATGEGCATCCGAAGACGTGGATGCGGGTGGTAACGGATCACATGGGTGITCTCTGAGA 
) GCCACHTGIGGTTTCTICIGITTICAGGCCCCAGAAGCCCATCCCGGGAAGGAAAATGCGC 
pane ( | CAATICCTGIGGCTITGGCCCTATCTGTCCTACGTGGAGGCTGTGCCCATCCGCAA 
GGTCCAGGATGACACCAAAACCCTCATCAAGACGATTGTICACCAGGATCAATGACATCTCACACACGGTAGGGAG 


GGACTGGGAGATGAGGTAGAACTGTGGCCATCCCTGGGGGAACCCCCAGAGGEC 


Fig. 2: NCBI BLAST data (Query: Reference Exon 2 sequence data, Subject: 
Aligned sequence data of the sample) 


AGTGCCTTTCATTACTGTCATTTCTAGACAATGAATTGTCTTTGAGGAGATGATAGCCAT 


PETED UTED UE EP EPP UP EP 
AGTGCCTTTCATTACTGTCATTTCTAGACAATGAATTGTCTTTIGAGGAGATGATAGCCAT 


GGCAGACAGCAAATCT 


PEVTET EET EPP ETE CETTE TEEPE EEE 
GGCAGACAGCAAATCTCATTGTTATCTGCATCCGAAGACGTGGATGCGGGTGGTAACGGA 


TETED UTEP EEE E TEE EE ETE EEE TEE 
TCACATGGGTGTTCTCTGAGATCGGCGACGGGCCACGTGTIGGTTTCTTCTGTTTTCAGGC 


AEC 
CATCCCGGGAAGGAARATGCGCTGTGGACCCCTIGTACCAATTCCTGTGGC 


PPTEV UT EU EEE EE UP ETE 
TTTGGCCCTATCTGTICCTACGTGGAGGCTGTGCCCATCCGCAAGGTCCAGGATGACACCA 


TAGGGAGGG 


PUPPET EET EE PEE ET PEPE EEE 
ARACCCTCATCAAGACGATTGTCACCAGGATCAATGACATCTCACACACGGTAGGGAGGG 


ACTGGGAGACGAGGTAGAACCGTGGCCATTICTGTGGGGGA--CCCCAGAGG 410 


PO ee oe ee eee 
ACTGGGAGATGAGGTAGAACTGTGGCCA-TCCCIGGGGGAACCCCCAGAGG 3 


Red: Leptin gene coding region 
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Jersey cattle were reported by Konfortov 
et al. (2006). Recent research in India with Bos 
grunniens revealed 4 SNPs in both exon 2 and 
3 as reported by Dongre et al. (2008). Among 
all these results, there was change of amino 
acids in protein product in some cases (i.e. non- 
synonymous mutation). Whereas, in few cases, 
no changes in amino acids were observed (i.e. 
synonymous mutation ) 

Jhala et al. (2009) have reported three point 
mutations at 42, 44 and 250 nucleotide positions 
in Mehsana buffalo. Orru et al. (2011) reported 
a total of 10 SNP: 2 in intronic region, 2 in exon 
2 and 6 in exon 3 out of which three SNP 
(C978T; C3100T; A3157G) have been reported 
to affect muscle fat in Simmental bulls. Scata et 
al. (2012) reported eight SNPs in 5' flanking and 
exon | of Bubaline leptin gene :A83G, A90G, 
A121G G256T, A283G, G9S59T, A1010C and 
G1254A. Tanpure et al. (2012) reported five 
single nucleotide polymorphic sites 98, 111, 172, 
209 and 266 in intron 1, which correspond to 
three haplotypes in Mehsana buffalo. 


Sequence Homology of buffalo leptin gene 

with other species 

The nucleotide sequences deduced for the 
respective exons of leptin gene in Murrah 
buffaloes were arranged to represent the coding 
region and were compared with other leptin 
sequences of the related species available in 
NCBI Gene Bank viz Bovine, caprine using 
NCBI web site Basic Local Alignment Search 
Tool (BLAST). There exists 97%, 97%, 99%, 
98% and 80% sequence homology with Bos 
taurus, Bos indicus, Ovis aries, Capra hircus 
and Homo sapiens respectively. The higher 
sequence similarity of buffaloes with sheep and 
goat than cattle might be due to lower query 
coverage (47%) in these two species. 


CONCLUSION 

Bovine leptin gene specific primers 
amplified the buffalo leptin gene and PCR 
amplification yielded an amplified product of 
exon 2 and exon 3 of bubaline leptin gene. Eight 
SNPs were identified while comparing bubaline 
leptin gene with that of other species, none of 
which resulted in change in amino acid sequence 
i.e. synonymous mutations. There exists 97%, 
97%, 99%, 98% and 80% sequence homology 
with Bos taurus, Bos indicus, Ovis aries, Capra 
hircus and Homo sapiens respectively. Because 
the leptin gene plays an important role in 
influencing economic traits in cattle, the novel 
detected single nucleotide polymorphisms might 
be used in association studies to assess their 
potential of being genetic markers for selection. 
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